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The JPL Low-Cost Silicon Solar Array Project is sponsored by the 
U.S. Department of Energy and forms part of the Solar Photovol- 
taic Conversion Program to initiate a major effort toward the 
development of low-cost solar arrays. This work was performed 
for the Jet Propulsion Laboratory, California Institute of 
Technology by agreement between NASA and DOE. 


This report was prepared as an account of work sponsored by the 
United States Government. Neither the United States nor the 
United States Department of Energy, nor any of their employees, 
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makes any warranty express or implied, or assumes any legal 
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SUMMARY 

The project began with techno-economic evaluations of 
MEPSDU processes, and then proceeded to a study of requirements 
for the longer range, potentially substantial improvement of 
the efficiency of Si solar cells. 

Of the MSP S DU processes investigated, ,tt was found that 
the Solarex metallization design and process selection should 
be modified to yield substantially higher output of the 10cm x 
10cm cells, while the Westinghouse design is extremely close 
to the optimum. In addition, further attention to the Solarex 
pn junction and base high/low junction formation processes 
could be beneficial. 

For the future efficiency improvement, it was found that 
refinement of the various minority carrier lifetime measurement 
methods is needed, as well as considerably increased sophisti- 
cation in the interpretation of the results of these methods . 

In addition, it was determined that further experimental inves- 
tigation of the Auger lifetime is needed, to conclusively deter 
mine the Auger coefficients for the direct Auger recombination 
at high majority carrier concentrations. This will determine 
the ultimately achievable efficiency of Si solar cells. If 
the Auger coefficients should be substantially lower than pres- 
ently thought, more attention may have to be given to bandgap 
narrowing. Finally, more needs to be known about the effects 
of various device processes on the minority carrier lifetime 
ultimately existing in the diverse layers of the device. 

Further, as one of the more frequently used methods for 
minority carrier lifetime measurement, the storage delay time 
(diode reverse recovery time) method was theoretically and 
experimentally investigated with respect to its limitations and 
its spurious effects which could lead to misinterpretation of 
the results. It was determined that two types of influences 
can alter the effects being measured and interpreted for the 
minority carrier lifetime of a given part of the sample: 
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influences from outside the part being measured (e.g. the base) , 
such as excess currents due to recombination in depletion regions, 
or currents due to changes of the depletion layer charge in 
transient methods, and secondly influences from the respective 
part itself, such as surface recombination of the back surface 
of this layer. A great deal needs therefore to be known of the 
structure and characteristics of the sample to be tested, before 
any method can be successfully applied and its results appro- 
priately interpreted. For the storage delay time method, these 
precautions have been delineated in more detail. 
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I . INTRODUCTION 


This project started under the LSA program and continued 
under the FSA project. Consequently, it was originally more 
closely coupled to the MEPSDU program, with technical and 
economic assessments of the various process approaches, and 
a particular view towards the avoidance of unnecessary solar 
cell performance losses in the process sequences applied. 

With the transition to the FSA program, the viewpoint 
shifted more towards future, larger efficiency improvements, 
and to the problems which stand in the way towards achieving 
these improvements. Consequently, attention was focused on 
the relationship between achievable minority carrier lifetimes 
and dopant concentration, with a particular view towards direct 
Auger recombination. As an outgrowth of this study, a closer 
look was taken at the problems and limitations of minority 
carrier lifetime measurement methods. 
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X. ASSESSMENT OF MEPSDU PROCESSES. 

A. Processes of the Solarex MEPSDU Sequence. 

The Solarex reports on the MEPSDU Contract DOE/JPL 955902 
give indications that a considerable improvement of the pro- 
cesses selected for the MEP.W0 sequence is possible. The aim 
of such improvement should be a substantial increase in cell 
efficiency, and consequently the module power output, while 
incurring not more than a moderate rise of the production 
costs. Thus, the cost-effectiveness of the process sequence 
could be significantly increased. 

Three areas have been noted to merit closer investigation 
in this respect: the contact design and metallization process, 

the pn (front) junction formation process, and the high low 
(back) junction formation process. 
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1* Front Metallization, 

In reviewing the Solarex MEPSDU cell design and process 
sequence# the item which seems to merit attention most is the 
front metallization system. According to the Solarex data^ 
which upon cursory check appear correct, 0.337W are lost due 
to front metal shading (0.139W grid line shading and 0.032W 
bus line shading) and to ohmic losses (0.068W in the diffused 
region, 0,014 in the grid lines, and 0.084W in the bus lines), 
while the cell output is only 0. 8-0. 9W, This amounts to a 
loss, just due to the front metal design, of 27 to 30% of the 
potentially available cell output if the shading and ohmic 
losses would be zero. 

Evidence of the insufficiency of the front contact pattern 
design and the metal selection is also contained in reference 2 , 
which describes several metallization experiments on 10 cm x 10 cm 
semi crystal line cells which were co-diffused in lot 10C with 
2 cm x 2 cm cells in a belt furnace. The small cells had TiPdAg 
metallization, the large cells the MEPSDU metallization process. 
Neither cells were AR coated. The small cells had an average 
efficiency of 6.7%, the large ones of 4.3%. Not all of this 
difference may be ascribed to the front metallization, since 
V QC and j sc also differed substantially. Also, as previously 
shown, using solder for the conduction layer has to be expen- 
sive, and is listed at $0.0216/W material content in reference 
3. In addition, this pricing is based on 15% efficient 
encapsulated cells , which is not likely to be attained with 
such large shading and ohmic losses, A finer-line grid, and a 
copper conduction layer will be needed to attain combined 
shading and ohmic losses as low as 12% on such large cells 
(10 cm x 10 cm) , and wire-built-up bus lines to reduce these 

losses further to about 5%. Of concern is also the relatively 

2 (4) 

high contact resistance of 0.26 to 0.41 fj*cm , which should 
find further attention. 

Additional questions on the quality of the metallization 

( 5 ) 

arise from the wide spread of the tab pull test data 


which ranged from a low of 100-200g (neglecting lower values 
from problem runs) Pj a high of about 1000 to 2000g, with 
standard deviations of *\bout 300 to 500g. It is not clear from 
the report whether the cells were solder dipped before the 
pull-tests which could provide some effect of sintering and 
strengthening the contact, but could also add the beam strength 
of the thick solder layer to the tab solder joint. The contact 
quality may also be questioned on the basis of the somewhat 
large values of contact resistance found, although these might 
be based on problems with the measurement method. 


2. PN-junction Formation Process. 

The second area where further investigations might be 
fruitful, is that of pn- junction formation. The data pre- 
sented^ 3 * * 6 ' 7 ^ do not conclusively indicate that 890°c is the 
optimum diffusion temperature. In fact, temperatures lower 
than 890°C seem to provide a higher solar cell power output. 

In addition, the application of a metallization system with 
too high a series resistance tends to let a diffused layer 
of lower sheet resistance appear more favorable, than a more 
optimized metallization would. Lower diffuse! layer sheet 
resistance values are obtained at higher diffusion temperatures . 

3. Base High/Low Junction Formation, 

The final question concerns the optimization of the high/ 

low Cback) junction formation process. 1 1 The data given show 
only averages, obtained on single lots, and do not permit an 
evalUAvi.Cvn of the statistical significance of the results. 

Including standard deviations and some repetition of the 
experiments helps to convey confidence in the process control 

achieved. 
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B« Processes of the Westinghouse MEPSDU Sequence, 

1. Metallization. 

As the Westinghouse fanned grid line pattern "appears so 
obviously wrong", with very dense grid-line placement in some 
areas , and very sparse placement in others -(Fig. l) , it was 
thought to be interesting to evaluate the total losses connected 
with this front layer metallization design. It was first at- 
tempted to make some quick, first order approximations, but it 
turned out that their errors were too large for any credible 
estimation. Closer approximations had, therefore, to be made 
for both the grid line losses and the front layer {diffused 
layer), losses, while computations of the total available power 
and of the grid line shading losses were quite simple. For 
the front layer losses, two computations were carried out, 
involving different simplified current flow patterns. Except 
for the irregular shaped area below the shortest grid line 
(located nearest the edge with the interconnect pads) , all 
equations have been derived. 

The dimensions of the pattern have been scaled off a 

( 9 ) 

graph for the 2.5 cm x 9.8 cm cell. For the sheet resis- 
tance of the diffused layer and for the dimensions of the 
grid lines, data were obtained from Westinghouse.^ 10 ^ The 
sheet resistance of the diffused layer is 60 ohms . The grid 
line width is 25 jjm, and their thickness is, on average, 8 ym. 
a. Description of the Evaluation Method 

The fanned grid line pattern as a whole is a periodic 
repetition of the pattern of a "unit cell". Each unit cell is 
symmetrical around its centerline which extends in the longi- 
tudinal direction of the unit cell. The pattern on each unit 
cell is the mirror image of the other, so that it is ade- 
quate to analyze a half-unit cell (Fig. 1) , The areas on the 
diffused layer from which current is collected by the individual 
grid lines, are of triangular shape. In addition, the triangu- 
lar collecting areas on. the two sides of each metallic grid line 
are of different shapes (Fig. 2) , It is therefore necessary 
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to evaluate the power loss separately for each triangular area 

of the diffused layer. It is simple to determine, for each grid 

line, the collecting area on the front layer. It is slightly 

more involved to determine the current distribution in the grid 

lines as function of distance from the interconnect attachment 

point. The current density entering the grid lines increases 

with distance from their attachment points up to a maximum 

near the far end of the grid lines, in contrast to the case 

of a parallel grid line pattern, where the entering current 

density is constant. The most difficult part is the determination 

of the power loss in the triangular (and one pentagon) shaped 

areas of the diffused layer. 

aa. Determination of the Available Power 

For the calculation of the available peak power, the 

approximate collecting areas on each side of the grid lines were 

calculated, and multiplied with the estimated maximum power 
. ~2 

current density j of 30 mA cm and maximum power voltage 

V of 0.5v. The areas A were determined from 
mp 

A = i b • h ; (.1) 


b being the length of the base of the triangle and h the length 
of its height. Thus, the available power is given by: 


avail 


h 




( 2 ) 


ab. Grid Line Shading 

The area shaded by a gridline is simply the product of the 
length of the grid line and its width T^. For the grid lines 
close to the centerline of the unit field, the angle between the 
centerline and the grid line is very small, so that the length 
of the grid line can be approximated by the longitudinal 
dimension W of the unit cell: 



W • T 


GL ; 


(3) 


ac. Determination of Joule Losses in the Grid Lines 

To determine the current distribution I (x) in the grid 
lines , a fundamental assumption had to be made concerning the 
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current distribution in the diffused layer. It was assumed 
that all current in the diffused layer flows normal to the grid 
lines. This is strictly true in parallel grid line systems, 
except near the juncture of the grid lines with a bus line. In 
the triangular system, a slight current spreading towards the 
interconnect attachment point should be expected. The effect 
will be very small in the narrow triangles, but may be more 
pronounced in the few triangles of greater height. 

The fanned grid line pattern contains 2 types of trian- 
gles (triangles a and b of Fig. 2) . In triangles of type a, 
some collecting area extends beyond the length of the grid line. 
For this part of the collecting area, it has been assumed that 
the current collected by it enters at the end of the grid line, 
and that this excess current compensates the reduced current 
entering the grid line from triangle b. 

For a reasonable approximation, it is thus assumed that 
both triangles are terminated by normal h to the grid line, 
which is located at its end. The current density entering 
the grid line at the distance x from the interconnect attachment 
area from one of these triangles with height h^ is then given by 



(4) 

(5) 


The power dissipated in a grid line element of length dx is 

given by: 


dP GL = I (x ; 2 dR ? 


where 


dR = R 


dx 

sh ' GL * T GL 5 * 7 


( 6 ) 

(7) 
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with R sh GL being the sheet resistance of the grid line metal- 
lization, and T gl the grid line width. Thus, the total power 
dissipated in a grid line is : 


P GL j mp~ R sh,GL 


(hi+h 2 ) 


GL 



which includes the contribution from both triangles of heights 
hj^ and h 2 , respectively. This yields 



-L i 2 r wL 

15 J mp K sh,GL T__ 

(a ij 


h l +h 2 


This compares with: 


(9) 


p as ± -i - c* . 

GL 7 J mp sh ,GL T gl & ' 

for the parallel grid line of uniform width, 
spacing S. 


( 10 ) 

with grid line 
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ad, Determination of Joule Losses in the Front Layer 

For the purpose of determining the Joule losses in the 
front (diffused) layer, two approximations to the actual 
value of this power loss were used. The more optimistic 
value resulting from those approximations has been used for 
the tabulation of the data. 

The analytical problem arises from the shape of the 
collecting area, which results in a case in which the equi- 
potential lines will generally not be parallel to each other 
or to the grid line (Fig, 3). Consequently, the current flow 
which occurs normal to the equi^potential lines, will not be 
normal to the grid line along its whole length, as it is in the 
parallel grid line case. Making the assumption that the 
current flow would be normal to the grid line (Fig. 4) , leads 
to a higher estimate of the power loss than will actually occur. 
In this triangular pattern, a current flow normal to the grid 
line would result in a very non-uniform current density along 
lines which are parallel to the grid lines. Consequently, the 
current will redistribute itself to paths which will be slightly 
longer, but which will result in a slightly more uniform current 
density (Fig. 3) . As this case is mathematically difficult to 
treat and not analytically solvable, the first estimation 
method assumes current paths normal to the grid lines. In 
contrast, the second method assumes that the current spreads 
out uniformly along lines parallel to the grid lines (Fig. 5) . 
The second method yields the lower estimates of the power loss. 

In the first approximation model, the current flow 
in a strip of width dx of the front layer (diffused layer) , 
oriented normal to the grid line and located at the distance 
x from the attachment point of the grid line, is given by: 

1 mp (x ' y) = j mp U “ y) dx '• 1 (11) 

The incremental power discipated in the length element 
dy located at y in this strip of front layer with the sheet 
resistance R g ^ is: 
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X (x,y) * dR PT ; 
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" ^mp ^ ~ dx ' R sh,FL * dx 
The total power dissipated in the strip is then : 

d W x) " j mp 2R sh,FL dx fu - y) 2 dy 


1 . 2„ .3 , 

" 3* D mp R sh,FL % dx 


As £ is a function of x: 


* ^ ~ W~x h ' 


and Ux) = — • h 
X 1 

x < W-x^ x > W-x.^ 


it is 


1.2 3 

T D mp R sh ,PL h 


(W-x x ) 


w 

| (W-x] 


P FL = h ^mpXh.FL h3 ”' (15 > 

The derivation which is general with respect to the shape 
of the triangle, and its results (eq. 15) show that the 
Joule losses are independent of the particular shape of the 
triangle . 

For the second type of approximation, a uniform current 
distribution is assumed to exist in a strip parallel to the 
grid line (Fig. 4) . Thus, the power dissipated in the strip 
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located at. the distance y from the grid line is expressed 
by: 


d^, L (y) - I(y) 2 dR(y) (16) 

As the current path is still assumed to be normal to the 
grid line, it is: 


and 


with: 


Thus : 


and: 


dR(y) 


m 


R sh,FL 



Ky) - j mp \ £ty) (h-y) ; 
£(y) «* W ; 


1 


dP Fl/ y) ~ 1 ^mp R sh,FL ^ y) * h ~ y) dy? 


p 2 r h 

*FL T J mp K sh,FL h 


/ 

o 


(h-y) J dy? 


P PL = ET Jmp 2R sh,F L • h " »» 


(17) 

(18) 


(19) 


( 20 ) 


The quality of the approximations is of some concern. 
For the parallel grid line pattern, the front layer Joule 
loss had been found to be: 


FL 


1 

= 1 


R 

mp sh,FL 


(I) 


W; 


( 21 ) 


for the rectangle of the front layer which is located on one 

S 

side of the grid line. The length ^ corresponds to the height 
h of the triangle on one side of the grid line in the fanned 

c 

pattern. As the total current from the rectangle is j W, 

, ^mp c 

but that from the triangle only j mp j W, a factor of ? for the 
square of the total current should be expected, and consequently, 
for the triangle, a Joule loss of: 
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P FL ** IT ^mp R sh,FL h W; 
as in equation (16) , 

For the derivation of eq. (16) , however, the assumptions 
of current flow normal to the grid line, and of a very non- 
uniform current distribution were made. It can be easily 
shown, that for a linear current distribution of the form: 

+ lixfl dx, 

■ 

< 22 > 

with the total current I (y) , the value of the square increases 
with the square of the slope a of the distribution. Thus, the 
square of the current is the smallest for the distribution 
assumed to be uniform according to Fig. 5. In this case, 
however, the current paths are longer than the assumed paths 
normal to the grid lines, so that eq, (20) would represent a 
lower limit for the actual Joule losses of the front layer. 
Equation (16) , on the other hand, will probably be an upper 
limit, as the current distribution will not be as non-uniform 
as assumed there, but some of this reduction of the losses 
will be compensated by the actual longer current flow paths. 

In any case, for the evaluation of the grid line pattern, 
the more optimistic values resulting from application of 
eg. (20) were used. 


I (x,y) dx » ax 


a: 
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b. Results 

The results of the calculations of the grid line losses 
for each grid line in the Westinghouse fanned grid line 
pattern, based on equations (2) , (3) , (9) , and (20) , are given 
in Table I, as well as their totals. It is seen that the 
approach of evaluating areas around all grid lines separately 
accounts for 18.08 mW of available power , while the avail- 
able power based on the whole half-unit field area should 
have been 18.75 mW. The computed total relative power loss 
of 4.7% is based on the available power value of 18,75 mW. 

For comparison purposes, a parallel grid line system 
has been evaluated with the same technology-based assump- 
tions: sheet resistances of the diffused layer of 60 ft, and 

of the grid lines (Cu, 8 ym thick) of 2.1 • 10~ 3 fl, and 
grid lines with uniform width of 25 ym. A four-line pattern 
on the half-unit field was found optimum under these con- 
straints, with the geometry and the losses shown in Table II, 
These losses include both the shading and the Joule losses 
of a bus line, which should give a more direct redundant 
connection to the adjacent unit fields, and consequently a 
connection of lower resistance, than the continuation of the 
grid lines into the adjacent unit fields does in the fanned 
pattern. The parallel grid pattern also includes an inter- 
connect attachment area identical to' that of the fanned 
pattern. 

In general, it can be concluded that the Westinghouse 
fanned grid line pattern design is a true low-loss metalliza- 
tion design. Nevertheless, a parallel grid line pattern can 
be used which performs at least equally well as the fanned 
pattern applying the same contact metallization technology 
constraints, Consequenlty , there is no performance uniqueness 
achieved with the fanned pattern. 
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4 Grid Lines 
25 p x 8 pirn 

1 Bus Line 

75 pm x 0.45 cm 

^ Contact Pad 

0.1 cm x 0.05 cm 


^ avail " 18 ‘ 75 mW 
Shading Losses: 

Grid Lines 
Bus Line 


0.3664 mW 
0.05 21 mW 



Contact Pad 

0.075 

mW 


Total : 

0.4935 

mW 

Joule 

Losses : 




Front Layer 

0.0763 

mW 


Grid Lines 

0.2186 

mW 


Bus Line 

0.0441 

mW 


Total : 

0.3395 

mW 

Total 

Losses : 

0.833 

mW 


2. 63 % 


1.81% 


4.45% 
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C . A Recommendation for Front Layer and Grid Lino Pattern 
Design ’ 1 ' 

An optimized grid line pattern design depends on the 
sheet resistance of the front layer, as it results from the 
diffusion or ion implantation process. Frequently, a partic- 
ular grid pattern design is chosen based on the thought that 
this pattern might be easily fabricated in an available pro- 
cess. Subsequently, the cell performance is optimized by an 
empirical process of varying time and temperature in the 
diffusion process . This approach generally does not lead to 
an optimum design, since a compromise is struck between the 
influence of the front layer design on collection efficiency, 
open circuit voltage and fill factor. However, only the 
latter is influenced by the series resistance which results 
from the sheet resistance of the front layer and from the grid 
pattern design. The proper approach would be to first opti- 
mize the front layer formation process to attain the highest 
product of collection efficiency and open circuit voltage. In 
a second step then, a grid line structure should be designed 
so as to minimize the series resistance losses from both the 
front layer and the grid line pattern, and simultaneously the 
shading losses from the grid pattern. 
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XI. SOLAR CELL DESIGN CONSIDERATIONS FOR FURTHER EFFICIENCY 

IMPROVEMENT 

A. General. 

All current cells employ, in the front layer, surface 
concentrations of the dopant impurity in the 10 cm range. 

It is well known that, in this concentration range, Auger 
recombination limits the minority carrier lifetime to values 
in the nanosecond range, while bandgap narrowing increases 
the saturation current and consequently reduces the open cir- 
cuit voltage. The first is the major reason why present 
diffused layers have to be made so extremely thin, the second 
results in the present limitation of the open circuit voltages 
to the low six-hundred millivolt range, and both contribute 
to the formation of the so-called "dead layer" at the surface 
of the front layer. 

A second part of the solar cell where heavy doping effects 
are frequently encountered, is in the so-called "BSF layer", 
if it is formed by diffusion or ion implatation. The BSF 
layers generated by these processes have been found to ba far 
less effective than those resulting from the aluminum alloying 
process, primarily because of the heavy doping associated 
with the common use of these two processes, although the gener- 
ally much smaller thickness of the highly doped layer between 
the high/low junction and the back contact which results from 
formation by these processes , would also reduce the effective- 
ness of the "BSF structure". However, with the limitation 
of the diffusion length to small values commensurate with the 
heavy doping, a thicker layer would not be more effective in 

this case. In the aluminum alloying process, the impurity con- 

18 — 3 

centration apparently is limited to 2*10 cm by the solid 
solubility of the aluminum in silicon, with a consequently 
longer diffusion length, to which the greater thickness resulting 
from the empirically found alloying process is better matched. 

For a future silicon solar cell of improved performance, 
heavy doping effects will need to be avoided in any part of 
the device. An exception might be a very thin layer just 
below the metallization, in order to facilitate the formation 
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of an ohmic contact. This very heavily doped and degenerate 
layer would then, in effect, form part of the ohmic contact, 
rather than an active semiconductor region. 

With the reduced doping in the front layer, a resistivity 
up to 2 orders of magnitude higher than currently experienced 
will have to be expected. One order of magnitude in sheet 
resistance should be regainable by increasing the thickness 
of the front layer to 1 to 3 pm, as the avoidance of the heavy 
doping effects should yield a minority carrier lifetime 2 to 4 
orders of magnitude higher than experienced in the present 
devices. This will still result in a front layer with a sheet 
resistance of 300-600 ohms. The proper reduction of the series 
resistance losses thus will require a fine line grid of about 
10-25 jjm line width with a line spacing of several hundred ym 
to a millimeter. The fabrication of such grid structures is 
well within the capabilities of today's technology, but may 
require the application of photolighography . 

In addition to avoidance of the heavy doping effects, 
there will be a need for control of the surface recombination 
velocity on the open front surface. While a number of methods 
for surface passivation have been explored and applied in 
industry, it is not clear what effective surface recombination 
velocities are actually obtained. A question also arises rela- 
tive to the results of some of the measurement methods applied. 

It will be necessary to pay more attention to the analysis and 
possibly the improvement of methods used for the determination 
of the front surface recombination velocity, which will be 
needed for the development of appropriate processes for its 
control . 

There is another large area of a dependence of the solar 
pell performance achieved on the processes applied. This con- 
cerns the influence of the high temperature processes on the 
material parameters actually attained. A part of this question 
again concerns the front layer,; this time with respect to the dif- 
fusion lengths which actually can be obtained by the processes 
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of diffusion, ion implantation, or CVD epitaxy. At least some 
of the processes used so far do not seem to be capable of 
attaining diffusion lengths comparable to those obtained in 
Czoehralski or float-zone grown material of the same resistivity. 
There is the additional question of the introduction of defect 
centers which act as recombination centers both in the deple- 
tion region and the base region near the depletion region. In 

the depletion region, such centers would cause an excess dark 

( 11 ) 

current, as described by Sah, Noyce, and Shockley . This 
excess dark current has found some study in the past and will 
require further investigation. If this dark current due to 
recombination in the depletion region is too large, it will 
degrade the fill factor. However, if the excess currents are 
negligibly small in the range of the IV-characteristic in which 
the current densities correspond to the maximum power in normal 
operating conditions, then a further reduction of these excess 
currents would not improve solar cell performance. 

Some evidence has also been reported in the past that 

defects resulting from diffusion are introduced into the base 

layer near the depletion region. Findings have recently been 

reported by NASA-Lewis, that the mobility in the base may be 

substantially lower near the depletion region than it is at a 

(i 2 ) 

greater distance from the depltion region . This would be 
a surprising effect, since the mobility is generally less sen- 
sitive to defects than the minority carrier lifetime is . 
Consequently, if deleterious effects on the mobility have been 
discovered, it should be assumed that the lifetime would be 
even more strongly effected, with an influence both on the 
collection efficiency and the open circuit voltage. The NASA- 
Lewis people properly state that a mobility degradation near 
the depletion region would improve the open circuit voltage, 
but an accompanying influence on the minority carrier lifetime 
should well counterbalance this effect. 


25 


pm? ^ 

OF .POOR QTO 



FIG. 6 


26 


MAJORITY CARRIER CONCENTRATION p(cm 


B. The Dependence of Minority Carrier Lifetime on Carrier 
Concentration. 

The values of the minority carrier lifetime, which can 
be expected to be obtained as function of the majority carrier 
concentration in silicon single crystal material were re- 
investigated in some detail through a literature study. Of 
particular interest are the limits of the achievable lifetime, 
including the question of direct band-to-band Auger recombina- 
tion. As one result of this study, the data from various 
authors of minority carrier lifetime measurements as function 
of majority carrier concentration in p-type silicon were 
plotted together on the same graph (Fig. 6). These data 
include those by lies and Socloff^^, by Kendall in their 

interpretation by Fossum^ 15 ! by Dziewior and Schmid \ by 

(17) i (18) 

Fischer and Pschunder , and by Mertens et. al. 

Particular attention was paid to the data which were 
thought to represent direct band-to-band Auger recombination. 

This direct recombination does not involve any traps or 
recombination centers , and is therefore presently considered 
as the physical mechanism which fundamentally limits the minority 
carrier lifetime at high majority carrier concentrations. Auger 
recombination principally involves two majority carriers and 
one minority carrier, and its rate therefore depends by the 
inverse square law on the majority carrier concentration. The 
effective minority carrier lifetime, t which is the result 
of two independent recombination mechanisms with the minority 
carrier lifetimes x x and 
by: 


T 2 occuring simultaneously, is given 


1 


eff 


t 


(23) 


2 


Auger recombination, with its 1/p dependence becomes the 
dominant recombination mechanism at adequately high majority 
carrier concentrations , and then determines the minority 


carrier lifetime x 


eff * 


As direct Auger recombination is a 
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fundamental mechanism not affected by defects or other process 
influences except the doping, it will determine the maximum 
attainable, or limiting minority carrier lifetimes. The Auger 
recombination is generally characterized by the Auger recom- 
bination coefficient C which is defined by the relationship 



(24) 


for p-type material. A value for this coefficient has been 
theoretically found by tedious quantum-mechanical con- 
siderations, involving a number of assumptions. This theo- 
retical value is therefore somewhat suspect. Dziewior and 

/t /r v 

Schmid. ' have experimentally investigated Auger recombina- 
tion in silicon by measuring the pulse of recombination 
radiation emanating from a wafer subsequent to the generation 
of minority carriers in it by the absorption of a strong photon 

pulse from a laser. Their lifetime data fall particularly well 

2 — 3 2 

onto a 1/p line, and yield an Auger coefficient 0^=9. 9 x 10 , 

in good agreement with the theoretically obtained value . A 

coefficient of this magnitude leads to the onset of the dominance 

of Auger recombination generally in the impurity concentration 

range of 10 17 to 10 18 cm“ 3 . 

Other workers, including lies and Soclof f ^^Fossum^ 15 ^ 

(19) 

and Possin et. al. , have been discussing relationships 
between minority carrier lifetime and majority carrier concen- 
tration, which follow a ^ dependency. While these authors do 
not try to relate the — dependency to any particular physical 

hr 

model of recombination, Figure 6 shows that these dependencies 

may primarily apply in the majority carrier concentration 

■16 18 “• 3 

range between 10 and 10 cm , that is below the range in 

which Auger recombination is dominant. A — dependency of 

P 

minority carrier lifetime can be explained by Schockley-Read- 
Hall recombination if the recombination center is a very 
shallow trap. Shallow traps, on the other hand, are usually 
not considered as being very effective recombination centers, 


so that this interpretation would require their presence in a 

high concentration. The levelling off of the minority carrier 

lifetimes towards still lower carrier concentrations, that is 
16 -3 

below 10 cm which is observed both in the Iles-Soclof f and 
the Kendall data, could indicate that there is also a deep 
recombination center present, which would result in a recombi- 
nation rate independent of the majority carrier concentration. 
This deep trap would have to be there in a lower concentration, 
so that it will determine the effective minority carrier life- 
time only in the higher resistivity samples. This "apparent 
saturation" of the lifetime in the higher resistivity samples, 
which is seen to occur at differing values of the lifetime, 
could be explained by differing concentrations of the deep level 
trap present in the various samples. 

Another interesting point concerns discrepencies in the 
experimental data at the high majority carrier concentrations. 

As explained above, Auger recombination should be the mechanism 
that poses the ultimate limitation to minority carrier lifetime 
at the high majority carrier concentrations. However, the data 

/TON 

from lies and Soclof f as well as those of Mertens v 

substantially exceed the Auger limit found by Dziewior and 

Schmid in their carefully executed experimentation. The Iles- 

Socloff data are quite consistent with those of Mertens , and 

yield minority carrier lifetimes a factor of 20 larger than 

those thought possible according to the Dziewior and Schmid data. 

In addition, an investigation of high cower transistor and 

( 19) 

thyristor structures by Possin, Adler, and Baliga has found 
that the exact modelling which includes band-gap narrowing and 
Auger recombination, cannot explain the experimentally obtained 
data from collection efficiency measurements by the EBIC method 
carried out with a range of electron beam energies , The inter- 
pretation of the data would require the assumption of Auger 
recombination at least an order of magnitude smaller than that 
found by Dziewior and Schmid. 
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The question arises where the differences in the life- 
time data from low resistivity material come from? One answer 
could be doubts in the accuracy of the minority carrier life- 
time measurements. Based on the general problems with minority 
carrier lifetime measurements, and the usually experienced 
discrepeneies , such doubts are justified. A recent study by 
this investigator has found that most of the minority carrier 
lifetime measurement methods are afflicted with side effects, 
which are usually neglected in the development of the theoreti- 
cal background of the method, and which usually are also not 
accounted for by the experimenter who applies the method. 

This entire area of minority carrier lifetime measurements is 
also one which will need considerable further attention. 

While, in principle, all experimental data are somewhat 
suspect, it is interesting to observe that the data with the 
longer minority carrier lifetimes are all based on measure- 
ments on completed devices , and may therefore include a greater 
potential for unaccounted side effects. In contrast, the 
Dziewior and Schmid data have been obtained on wafers, without 
further processing or introduction of a potential barrier after 
the growth and the cutting of the crystal. Also, the investi- 
gation of the decay of the recombination radiation as performed 
by Dziewior and Schmid is thought to be less susceptible to side 
effects, than, e.g., the diode recovery method used by lies 
and Socloff. Also, the fact that the experimental results of 

Dziewior and Schmid show only very small deviations from the 
2 

1/p relationship conveys some confidence in the acauracy of these 
data. On the other hand, the three investigators using devices 
have applied completely different measurement methods as well as 
different device structures, and have still arrived at results in 
close agreement with the conclusion that the Auger coefficient 
should be more than an order of magnitude smaller than that found 
by Dziewior and Schmid. This agreement does not seem to permit 
neglecting these data offhand. In fact, Possin et. al. < ’ 20 ^ 
outline a possible interpretation for the observed discrepancy: 
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It is known that trap Assisted Auger recombination also exists, 
besides the direct band-to-band recombination discussed above. 
Both types of Auger recombination show the 1/p dependency. Thus, 
Dziewior and Schmid may have had samples in which trap assisted 
recombination was the dominant effect, rather than the direct 
recombination which they thought they were observing. In the 
preparation of devices, significant gettering can take place, 
so that the trap concentrations in the device samples could "have 
been adequately suppressed so that the authors who investigated 
the devices and observed the longer minority carrier lifetimes, 
may actually have seen the direct band-to-band Auger recombina- 
tion. Clearly, both explanations are sheer speculation at 
this time. 


III. EVALUATION OP MINORITY CARRIER LIFETIME MEASUREMENT 

m wmur. 


A. General 

The preceding section has directed attention to the 
limitations of minority carrier lifetime measurements# and to 
the interpretation of data obtained in such measurements# in 
connection with the literature data on Auger recombination. 

There is an additional need for reliable minority carrier life- 
time measurements, namely for the determination of the minority 
carrier lifetimes actually achieved in the various layers of 
high efficiency solar cells# if the processing of these cells 
is to be controlled so as to yield the cell performance which 
is being predicted by modelling. 

A large number of minority carrier lifetime or diffusion 
length measurement methods is in existence . Some of these 
methods are only minor variations from another method. In 
principle# the measurement methods can be divided into several 
categories. The first division may be into methods which are 
based on steady state (do) phenomena, those, which are based 
on periodic (ac) phenomena# and those which are based on tran- 
sient effects. The steady state methods generally determine 
the diffusion length rather than the minority carrier lifetime. 
The transient phenomena include effects such as recovery after 
a time period of excess minority carrier injection. A second 
broad categorization is based on the type of sample required. 

The first category uses a piece of plain semiconducting material, 
generally with at least 2 ohmic contacts. An example of such 
a method is the photo- conductive decay method. The second group 
requires the existence of a potential barrier in the sample, 
such as a pn junction or a Schottky barrier. 

All the methods have been originally developed for a 
simplified case, which is relatively easy to model, and where 
the mathematics is tractable. The assumptions made to facili- 
tate such modelling vary somewhat from method to method, as do 
the effects which have been neglected.. The problem is# that the 
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methods are frequently rather blindly applied, without attention 
to their limitations, or the degree of conformity to the assump- 
tions in the particular sample being investigated. The literature 
is full of examples where such neglect is evident. 

An effort has therefore been started to determine the 
limitations of various minority carrier lifetime measurement 
methods, and of the types and magnitudes of corrections to be 
applied, as far as such corrections have been derived, or can 
be derived. It can readily be seen that it is necessary for 
the investigator of the minority carrier lifetime in a particu- 
lar sample, to know the structure and the characteristics of 
the sample in great detail before he can properly interpret the 
data which result from the application of a particular minority 
carrier lifetime measurement method. Only with such knowledge 
can he apply the method in a range where disturbing effects 
are least influential, and where he can interpret the results 
of his measurements most meaningfully. Because of this consid- 
erable influence of extraneous effects, at least 2 completely 
different methods should be applied to any particular sample, 
and agreement between the minority carrier lifetime data re- 
sulting from these two different methods will be needed to 
attain credibility for the data. 

Of the numerous minority carrier lifetime measuring methods , 
the storage delay time (diode reverse recovery time) method has 
been studied in more detail for several reasons. It is a method 
which has been very commonly used because it is easy to apply 
with relatively inexpensive equipment. Many data in the lit- 
erature, including those in several of the "Auger lifetime" 
articles, were obtained by use of this method. Also, the equip- 
ment for this method is available in this investigator's labo- 
ratory, so that experiments can readily be carried out. The 
method has, however, lately been touted as "unreliable", 
although its limitations and problems of application and inter- 
pretation may not differ much from those of most other methods. 


Further, a Master degree student of this Principal 

investigator has last year expanded the theory for the diode 

reverse recovery time or storage delay time (SDT) method for 

minority carrier lifetime determination to the case of arbi"* 

trary base widths with arbitrary surface recombination velocity 

at the plane boundaring the base opposite the pn- junction. 

All previous publications have dealt only with the cases of 

(21 22 ) 

infinitely thick base regions ' , or of infinite surface 

( 23 - 26 ) 

recombination velocity with arbitrary base widths 

However, the modern solar cells, containing in the base 
region a high- low junction and a third, more heavily doped layer 
of varying width, can exhibit a wide range of transport veloc- 
ities through the interface plane between the first base layer 
and the high-low junction, depending on the design and quality 
of the high-low junction and the third layer. This transport 
velocity is equivalent to a surface recombination velocity at 
the back surface of a single layer base region. 

In the real cases, both the bulk minority carrier lifetime 
and the transport velocity at the rear of the base layer are 
unknown. It will then be necessary to prepare several (at least 
two) samples with differing thickness of the base layer, but 
otherwise identical properties, to determine both the lifetime 
and the transport velocity. 

With this new capability and the other positive attributes 


of the SDT method, a closer look seems Indicated to determine 
its limitations, its range of applicability, and the precautions 
to be taken in its use. 
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B. The Storage Delay Time (Reverse Recovery Time) Method of 
Minority Carrier Lifetime Measurement. 

This method belongs to the group which utilizes a transient 
effect, and which reguires a sample with a potential barrier, 
usually a pn junction. A closely related method is the open 
circuit voltage decay (OCVD) method. In one variant of the 
latter, the same equipment is used as for the storage delay time 
(SDT) method, so that both methods can easily be used in con- 
junction with each other. Nevertheless, the OCVD method will be 
discussed separately. 

The SDT method is rather easy to use and not time consuming. 
It requires only a fast switching device, two power supplies, 
and an oscilloscope (Pig. 7) . For the fast switching device, 
a mercury wetted contact relay is frequently used, as it has no 
contact bounce. In another variant which is better suited for 
the measurement of very short storage delay times (< lys) , a 
power supply and a high power pulse generator are used. 

The experimental procedure involves the measurement of the 
storage delay time t on the oscilloscope display, after the 
forward bias injection current I^ and the reverse bias recovery 
current I have been appropriately adjusted (Pig. 8) . If dis- 
turbing electronic effects are present, such as significant 
overshoots, pulse rounding, or ringing, they are displayed and 
can either be abstracted by visual interpretation, if modest, or 
be eliminated by improved circuit design. 

The method has originally been derived ^ 21 ' 22 ^ under the 
assumption of an abrupt, highly unsymmetrical junction, which 
means that there is only insignificant injection into the more 
heavily doped side of the pn junction, and consequently no sig- 
nificant charge storage there. It has also been assumed that 
-the base region is infinitely thick, so that there is no influ- 
ence of the boundary surface to the base opposite the pn junction. 
As mentioned inthe introduction, a number of workers ^ 22 ~ 2 ^ have 
tried to expand the theoretical background to determine the re- 
lationship between the measured storage delay time and the 
minority carrier lifetime in the sample for the case of arbitrary 
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base thickness, but with infinite surface recombination velocity 
at the back surface, 

In the SDT method, a steady state excess minority carrier 
distribution is built up in the base region by injection of a 
current across the pn junction with forward bias applied. For 
the infinitely thick base, this distribution is a negative 
exponential, with the excess minority carrier concentration 
n p( x j afc i nterface between the depletion layer and the 
base (Xj determined by the forward bias voltage V existing 
across the depletion layer: 


S 


n (x . ) - n 

p 3 ,B po 


* n (e 

DO 


- 1 ) 


(25) 


After switching at time t=0 from the injection condition 
to the reverse bias power supply, the minority carrier distri- 
bution in the base is initially the same as before, except for 
an instantaneous change of the slope of the distribution at 
x=x. „ from negative to positive, and a consequent minute change 

J a . 

in n (x. _) . With n (x. n , t=0 ) « n <x. D t<0) immediately 

P, 3 / B PDf B / 

after switching (t= 0 + ) , the voltage across the depletion region 
has also still to be the same with respect to sign, and only 
minutely decreased in magnitude. 

As time progresses, the excess minority carrier distribution 
in the base will decrease for 3 reasons: 1. recombination in 

the bulk of the base; 2 . outflow of minority carriers into and 
across the depletion region beaause of the large electric field 
in this region (analogue: light generated current from base) ; 

and 3. outflow of minority carriers through the rear boundary of 
the base region, if the latter is less than a few diffusion 
lengths thick. Consequently, n^ (.x^ t) will decrease with 
time, and so will the bias voltage V. The outflow of carriers 
across the depletion region, and consequently the current I , is 
limited by the resistance R. As the forward bias across the 
depletion region is in the same direction as the power supply 
voltage V , the voltage drop across the resistance R will be: 
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Eventually, n (x . t ) = n will be reached, and 
P 3 t " s P° 

simultaneously V=0 . From this point in time (t « t g ) on, the 
junction goes into reverse bias, although there is still stored 
charge in the base (Fig. 9). Thus: 

V -V 

I (t > t ) = - < I CO < t < tj 

I decreases until it asymptotically reaches I Q . 

During the entire process, the distribution of the excess 
minority carriers in the base changes, depending on the rela- 
tive magnitudes of bulk recombination and outflows across the 
two base region boundary planes. In consideration of this 
change, the relationship between the storage delay time t and 

9 

the bulk minority carrier lifetime x n in the infinitely thick 
base has been found to be^^' . 

— y ' ts (27) 

1 + j 

(27) 

Streetman has shown that the mathematical derivation 
can be greatly simplified when the assumption is made, in the 
case of the infinitely thick base, that the distribution will not 
change in time, but rather stay exponential. Then: 

t = — = — r — t ( 28 ) 

n »„(i + !£) s 

x r 

For the usually applied case of I =1^, the factor on t is 4.5 

IT L b 

in eq. (.27), and 1.44 in eq. (28), which is a- very significant 
difference . 

The method always involves the measurement of the time 
from the moment of switching to the end of the constant current 
phase, that is the "storage delay time". For a given sample, the 



39 


ORIGINAL PAOS W 
OF POOR QUALITY 



Figure 9 


storage delay time varies according to the ratio of forward 
injection current to reverse recovery current. The relation- 
ship between the storage delay time and the ratio of the currents 
is given by eq. (27) only, if the diode is infinitely thick, 
and no other extraneous effects are present. 

The limitations of the method are manyfold. The first 

limitation to be eliminated had been the restriction to large 

base widths, albeit only for the case of infinite surface 

recombination velocity at the back surface. In the last year, 

as mentioned before, the theory was expanded to include arbitrary 

( 28 ) 

surface recombination velocity with arbitrary base widths 
(Fig. 10, 11). A problem is that, normally, the surface recombi- 
nation velocity at the back boundary of the base layer is just 
as little known as the minority carrier lifetime. However, by 
preparing samples of varying base thickness, e.g. on the same 
wafer and with the same process, it should be possible to 
determine both the surface recombination velocity and the dif- 
fusion length in the base layer simultaneously. This would be of 
considerable benefit for the evaluation of solar cells which in- 
clude high-low junctions. In these cells, it is important to 
determine the effectiveness of the high-low junction structure 
prepared, as well as the minority carrier lifetime effective in 
the base after all the applied processing. 

There are several effects which can influence the deter- 
mination of the minority carrier lifetime by the storage delay 
time (SDT) method. These effects can be grouped into two catego- 
ries: first are influences on what appears to be the storage 

delay time resulting from the charge stored in the base. These 
influences result from the opposite side of the depletion region, 
(as seen from the base) , from excess current due to recombina- 
tion in the bulk or on the surface of the depletion layer of the 
pn junction, and from currents which result from the change of 
the depletion layer charge with varying junction voltage. The 
second group of effects includes those which are connected with 
the base region itself, but give rise to a base region "relaxa- 
tion time" which does not equal the minority carrier lifetime 


41 


PIG. IQ Dependence of the ratio of the recovery time to the lifetime 
On the ratio Of the base width to the diffusion length for 
various values of the ratio of the surface recombination 
velocity to the bulk recombination velocity (the ratio of the 
diffusion constant to the diffusion length) for Jf s J r . 

(The data for Fig. 10 and 11 were recalculated, with higher 
accuracy, after completion of the contract, and new figures 
representing the new data, were inserted in this version of 
the report.) 









’-n the bulk of the base region. These influences are the trans- 
port velocity out of the rear boundary surface of the base layer 
(frequently a surface recombination 'Velocity) , and the surface 
recombination velocity at the circumferential surfaces of the 
base layer. 

The first severe criticism of the method was raised by 
( 29 ) 

Neugroschel et. al. who pointed out that any significant 
charge storage on the opposite side of the pn junction from the 
base (they call it "emitter" although it may not have emitter 
properties) will influence the relationship between the minority 
carrier lifetime in the base and the measured storage delay time. 
They point out that the method really involves the determination 
of a "relaxation time" which includes the charge decay both in 
the base and the "emitter", but omit mentioning any back surface 
effects. They properly point out that the derivation of the method 
includes only the diffusion component of the diode forward 
current, and that therefore the measurement should be performed 
only in the region where the diffusion component determines the 
IV characteristics of the device. In real diodes, excess current 
results from recombination at the surface and in the bulk of 
space charge layers, particularly the depletion layer of the pn 
junction, and it is usually the dominant current in the low- 
voltage regime of the forward IV characteristic. 

For a successful application of the method f therefore, the 

current voltage characteristic of the device to be tested has 

( 29 ) 

to be known and considered. Neugroschel et. al. prefer the 

OCVD method over the SDT method, since they feel that the former 
is self-indicating if it is not carried out in the regime where 
the diffusion component dominates. However, if the SDT method 
is carried out at current levels where the excess current is 
negligible, it is fully valid, despite a statement to the con- 
trary by Neugroschel et. al . Since the method is carried out 
with constant current flow, these current levels for 1^ and I 
can, for many diodes and for most solar cells, easily be chosen 
so that they fall more than one order of magnitude above the 
excess current, while remaining in the "low-level" injection regime. 
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In the forward bias injection phase, any recombination 
in the bulk or at the surface of the depletion layer results 
in an "excess current" eJ£C , by which the charge injected and 
stored in the base region is reduced from that expected to 
result from the measured forward injection current 

f ,meas 

In the reverse bias recovery phase, the highest value of the 
excess current I (t) will occur at the beginning of the 
storage delay phase, as the excess current increases with for- 
ward bias voltage, albeit less than the diffusion current does. 
If this amount of excess current is neglible, then the sub- 
sequently decreasing amounts will be negligible. It should be 

noted that, during this phase, the extraction current is in- 

£ 

dependent of the IV-characteristic, based solely on the stored 
charge and its distribution, but limited by the external re- 
sistor. Because of the still existing high carrier concentra- 
tion in the depletion region (forward bias!), there occurs 
still enhanced recombination there, which leads to a current 
in the direction opposite to that of the extraction current, 
and thus to a decreased value of I r « Consequently, the measured 
current values are; 

■^meas ^f * X f,exc ' 

and; 

r,meas “ r r,exe v ' ' 

As one adjusts I to be ecrual to I- , then 

J r,meas - f,meas 

obtain : 


( 29 ) 

( 30 ) 

one will 


I JT + I . + I „ ^(t) 

f 0X0 f £ 0X0 f JT 


> I- + X JZ t 

^ f exc,f ' 


so that too short a storage delay time will be obtained. If 
I f cannot be chosen large enough to make I F « I- and 
consequently negligible, knowledge of the IV characteristic 
can permit a correction for I , F by setting; 

0X0 / i ” 


45 


^r/meas * ^f^eas " ^exc^f * *£ ^ 

This may permit acceptable measurements to be made in the range 

where I - < I-, but not I * « I-* 
exc,f f f exc/f f 

A second assumption usually made is that the charge change 
in the depletion layer is small compared to the charge stored 
in the base. 

As the depletion layer charge increases only with the 
square root of the applied bias, its effect can be made negli- 
gible by going to high enough current levels for the storage 
delay time determination, just as was found to be the case with 
the excess currents due to recombination in the space charge 
layer. The magnitude of the depletion layer charge can be deter- 
mined either by a direct measurement of the junction capacitance, 
by calculation of the junction capacitance from known design 
parameters of the device - in either case at the values of 
applied bias at the beginning of the recovery phase and at its 

end - or from the decay phase following the recovery phase in 

(31) 

the SDT method as suggested by Kuno's relationships. 

Practical application of the SDT method should consequently 
be preceded by an analysis of the IV characteristic of the device 
to be tested (Pig. 12 as an example) and a determination of the 
depletion layer charge. Figure 13 shows four storage delay traces 
of the 1cm x 2 cm Si n/p solar cell whose IV characteristic is 
shown in Fig. 12, taken at four different injection levels 1^ meas 
for which the excess current is insignificant, as comparison with 
Fig. 12 indicates. At the three higher levels of the injection 
current, the same value of storage delay time is obtained, al- 
though the storage delay pulse becomes successively more rounded 
as the injection level is decreased. It may be noted that in 
all these cases, the low level injection condition is maintained. 
At the lowest injection level, however, the storage delay pulse 
is nearly completely obliterated by a "decay curve". This 
successive rounding and final obliteration is thought to result 
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from the depletion layer charge, although this point is not 
proven. It cannot be the result of the excess current, but 
the change in depletion layer charge is of the same magnitude 
as the charge stored in the base at the lower injection level. 

Finally, if the surface recombination velocity at the 
open surface at the circumference of the base layer is high, 
the excess minority carriers within approximately 1 diffusion 
length from this surface will diffuse to this sink rather than 
to the pn junction or the back surface. Thus, the existence 
of a high surface recombination velocity at the circumference 
of the base layer can have a significant influence on the 
relaxation time of the base layer, and consequently on the 
supposedly measured carrier lifetime for this layer. 
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1. Kuno's Modification of the Storage Delay Time Method. 

Kuno has tried to simplify the theory of the SDT method, 

(31) 

and to expand its application. His extensive use of the 
charge control approach has made the analytical treatment 
appear very elegant, but has led to oversimplifications which, 
in the view of this investigator, have made his particular 
solutions useless for minority carrier lifetime determination. 

Kuno has introduced the relaxation time x,- which takes into 

r 

account minority carrier bulk recombination in the base layer 

as well as transport out of this layer through its back surface. 

This relaxation time x^ is a constant during the steady state 

part of the injection phase, but will generally not be constant 

during the recovery phase, as the carrier distribution in the 

base layer changes during the recovery phase, so that the ratio 

of the average minority carrier concentration in the base layer 

to the concentration at its tear boundary surface will vary in 

time. Consequently, the ratio of the bulk recombination to 

the transport out of this surface will vary during this phase, 

and with it the relaxation time x^. The exception is, when 

only bulk recombination is present, so that the relaxation time 

is constant. Treating this relaxation time as a constant, Kuno 

obtained an analytical solution to a differential equation for 

the total stored charge Q(t) during the recovery phase with 

constant extraction current -I : 

r 


dQ Q_(t) 
dt x^ x r 

The solution is : 


(32) 


Q ( t) = (I f + I r ) x f exp(-t/x f ) - I fc x f (33). 

With a non-constant relaxation time, the time dependence of the 
total stored charge should differ from the exponential, with 
an analytical solution not possible. 
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The second problem arises from Kuno's definition of the end 

of the storage delay phase. The time t is generally defined 

s 

as the time from switching to reverse bias until the instant when 

the excess minority carrier density at the boundary surface 

between the base and the depletion region reaches zero. During 

this entire period t , the diffusion current due to the gradient 

s 

of the excess minority carrier density in the base sustains the 
current -I r , which is limited to this constant value by the 
external resistor R. Kuno defines the total charge stored in 
the base at time t by the expression: 

Q(t s > = I r T r (34) 

and replaces the proper boundary conditions 


p(Xj,t g ) = o ; 


(35a) 






(35b) 


by eq. (3 4) . He thus replaces the unknowns Q(t ) and t by 

the unknown time constant x . By this approach he arrives at 

a relationship between t and x f : 

s x 


t 


s 


Jin (1 + 


I f T r 

^) - £n(l + ~) 
r f 


(36) 


Plotting t against £n(l + •=■— ) should thus yield a straight line, 
® If 

with the slope x f and the intersect with the abscissa at •=— = 

, Tr 1 - L E' 

~r-, as seen in Fig. 14. 

T f 


The storage delay time was measured on a number of diodes 

I f • 

and solar cells with varying values of I f and =— . Plotting t 
If r 1 r s 

versus £n(l + j — ) has actually led to straight lines in many, 
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FIG. 14 



but not all cases, but the values of x r obtained have all been 
small compared to t^, or actually zero. Consequently, the 
second term in eq. (36) has always been essentially zero, re- 
ducing the relationship between and t g to the logarithmic 
relationship obtained by Streetman by assuming that the excess 
minority carrier distribution does not change during the recovery 
phase (eq. 28) . Consequently, the minority carrier lifetime 
determined thusly is approximately a factor of 3 too small com- 
pared to the one obtained under consideration of the changing 
excess minority carrier distribution. 

The problem evidently arises from Kuno's use of eq. (35) 

with various values of I , considering t as a constant. In 

r r 

reality, Q(t ) will not be proportional to I at varying values 

of I r , for values of t as usually defined, and as measured in 

the experiment. With T r = constant, a different t Kuno will 

be defined, which determines a Q(t v ,,„„) according to eq. (35), 

s f j\uno 

but this t g Kuno will not be measureable in the experiment. 

Kuno extended his theory to also include the decay phase 
which follows the storage delay phase. While neglecting the 
effects of the depletion layer capacitance Cj during the storage 
delay phase, as have the previous workers, he has included its 
effect during the decay phase. Thus, he arrives at the relation- 
ship : 


2.3 


T , RC . 
r + n 

1 + T /T* 

r f 


(37) 


for the measurable time t f which he defines as the interval from 

the end of the storage delay phase (t = t ) to the point where 

the current has been reduced to 0.1 • I . The capacitances 

obtained by this method were found to be about 50% higher than 

the zero bias capacitances calculated for the devices, while 
the actual effective capacitance during the decay phase should 
be less than the zero bias capacitance. 
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C. EXPERIMENTAL 

The equipment for the determination of the dark IV charac- 
teristic and that for the observation and measurement of the 
storage delay effect, as well as for the OCVD method, have been 
assembled, carefully checked out, and debugged. They have been 
applied to a number of diodes which show a significant storage 
delay effect, as well as to several two ohm centimeter n on p 
silicon solar cells which were left over from an earlier ex- 
periment, At least one of these solar cells had been irradiated 
with low energy protons. 

It is interesting to note that those of the diodes which 
exhibited significant storage delay effects, and which were 
available from supply, displayed strange IV characteristics. 
These characteristics seemed indicative of high level injection 
effects, as well as of narrow base effects. Apparently in con- 
sequence of the latter, the recovery characteristic shown in 
Fig. 15 was obtained, which exhibits an ideal recovery pulse 
with a negligible decay phase. 

The solar cells, all of approximately 1cm x 2cm junction 
area, of either planar or conventional construction, with 2 Qcm 
base layer and a p+ layer, when measured with I f - I = 250 or 
100 mA, showed a less square recovery pulse, and a significant 
decay phase (Fig. 13). At the current levels applied, both 
the excess currents and the depletion layer charge are an order 
of magnitude smaller than the injection current and the charge 
stored in the base, respectively (Fig. 12) . The "Kuno interpre- 
tation" of a number of storage delay time measurements on this 
cell at three levels of 1^, and with varying I^/I^ ratios, gave 
neither credible straight lines, nor the same slope or time 
constant T r for the three injection levels (Fig. 16) . In 
addition, the fall time t^ increased from 7 ys to 20 ys with 
increasing I^/I.^. ratio from 1 to 50, at 1^ = 250 mA, with a 
similar behavior observed at 1^ = 100 mA. 
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The smallest value of fall time t f in the decay phase , 

7 ys, together with the largest value of t , 1.8 ys, resulted 
in the capacitance value C.. = 44 nF. This value compares with 
the measured zero bias capacitance value of 68 nF, and the 
calculated value of 52 nF, for a rather good agreement in this 
particular case. 

The minority carrier lifetime in the base would be 13.5 ys 

according to the Kingston theory for the infinitely thick base. 

The Kuno interpretation gives <6.4 ys , and OCVD resulted in 

2 -T -l 

1.6 ys . With y n - 1200 cm V s , the diffusion length with 

13.5 ys lifetime would be L - 200 ym. With a cell thickness 

n 

of about 300 ym, and a high/low junction present, the ''infinite 
base" approximation is questionable. Such questioning is par- 
ticularly justified, as the measured saturation current density 
j Q ^ would demand a twofold larger diffusion length, or a five 
times larger minority carrier lifetime than obtained from the 
Kingston method. This difference is not readily explained. 
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XV. CONCLUSIONS 


A review of the MEPSDU metallization processes has 
indicated that the Solarex grid pattern design and process 
selection, in combination, lead to excessive Joule and shading 
losses, while the Westinghouse metallization system is as 
close to optimized as can be expected. In addition, it is 
not clear that the Solarex front layer diffusion process and 
the BSF formation process are actually optimized. 

A study of the minority carrier lifetimes, achieved in 
various samples, in dependence on the impurity concentration 
has led to three primary conclusions. First, further investi- 
gation and refinement of the various minority carrier lifetime 
measurement methods and of their interpretation is required. 
Second, a further study of Auger recombination is needed, until 
the Auger coefficients for both trap-assisted and direct re- 
combination are conclusively determined by application of 
differing methods. Third, some device processing methods are 
deleterious to the minority carrier lifetime, while others may 
include adequate amounts of getter ing so as to improve the 
minority carrier lifetime by an order of magnitude or more 
relative to that found in the as-grown and cut samples. Since 
the minority carrier lifetime is one of the key parameters which 
determine the efficiency available in silicon solar cells, more 
studies will have to be performed to understand the influences 
on the minority carrier lifetime which is exerted by the various 
processing methods available. 

An additional, but probably less urgent conclusion is, that 
the ultimate determination of the direct Auger recombination co- 
efficients will lead to a review of the optimum design for the 
solar cell. If the Auger coefficients should actually turn out 
to be more than an order of magnitude smaller than the presently 
accepted value, the design of the solar cells would be changed 
in the direction of utilizing higher impurity concentrations, 


r 




with corresponding adjustments in the device structure so as 
to achieve an additional efficiency gain. 

An investigation of minority carrier lifetime measure- 
ment methods has led to the conclusion that all methods are 
afflicted by 2 sources of errors: effects which lie outside 

of the portion (such as the base layer of a diode) of the 
sample for which one wants to determine the minority carrier 
lifetime; and effects which originate from the desired por- 
tion of the sample, but which alter the "signal" which would 
be attributable only to the bulk minority carrier lifetime. 

Most prevalent among the latter are surface recombination 
effects, while for measurements on devices with a potential 
barrier, excess currents and depletion layer charges are fre- 
quent sources of error from the first group. To arrive at 
credible data, a great deal needs to be known about the attributes 
of the sample to be tested, and considerable sophistication 
exercised in the execution and interpretation of the measure- 
ments . 
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